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X) . ABSTRACT 

<L> ■ 

pL| ' Aims. We investigate the detection of Cool Cores (CCs) in the distant galaxy cluster population, with the purpose of measuring the 

CC fraction out to redshift 0.7 < z < 1.4. Using a sample of nearby clusters spanning a wide range of morphologies, we define criteria 
to characterize cool cores, which are applicable to the high redshift sample. 

Methods. We analyzed azimuthally averaged surface brightness (SB) profiles using the known scaling relations and fitted sin- 
gle/double p models to the data. Additionally, we measured a surface brightness concentration, csb, as the ratio of the peak over 
the ambient SB. To verify that this is an unbiased parameter as a function of redshift, we developed a model independent "cloning" 
technique to simulate the nearby clusters as they would appear at the same redshifts and luminosities as those in the distant sample. 
This method is based on the application of the cosmological surface brightness dimming to high resolution Chandra images and as- 
sumes no intrinsic cluster evolution. A more physical parameterization of CC presence is obtained by computing the cooling time at 
a radius of 20 kpc from the cluster center. 
Results. The distribution of the SB concentration and the stacked radial profiles of the low-z sample, combined with published infor- 
mation on the CC properties of these clusters, show 3 degrees of SB cuspiness: non-CC, moderate and strong CC. The same analysis 
applied to the high-z clusters reveals two regimes: non-CC and moderate CC. The cooling time distribution corroborates this result by 
showing a strong negative correlation with c$b- 

Conclusions. Our analysis indicates a significant fraction of distant clusters harboring a moderate CC out to z= 1 .4, similar to those 
found in the local sample. The absence of strong cooling which we report is likely linked with a higher merger rate expected at redshift 
z > 0.7, and should also be related with the shorter age of distant clusters, implying less time to develop a cool core. 
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<0 1 . Introduction I2005I) . and the mass deposition rates (M) are much smaller than 

, expected (Edge & Frayer 2003). 

, Galaxy clusters are filled with hot, diffuse gas - the intra clus- 
ter medium (ICM) - which strongly emits in X-rays. While the 



The existence of a feedback mechanism which counteracts 



£*■ ■ , . , , j* . •! .• ■ ,, , j, .. cooling is now widely accepted and currently, the most debated 

.^h cluster dark matter distribution is generally described by its grav- . & . , . , J . F , . , . . 

S<! . itational potential (Navarro, Frenk & White 1997), the physical P 1 ™. 18 bBSb ^J c ^f gal f tlC nU ? le ! (A ?, N) (se f , reVleW 

U ' processes which govern the behavior of the ICM are complex ^ Fablan et ■MESS. Although several plausible models exist 

rn ' a • i a Z •<- <-: i u which may explain how this mechanism works (Churazov et al. 

L « , and include non-gravitational phenomena. Cool cores are a man- — — . / F n . , , ., 

ifestation of these intricate processes. 2000 ' Bm 8f n & , Kalser 2002 ^ and there 1 ls evldence for enou 8 h 

. j. . , t . , T „ t, i lr , 0/l | t, , • output mechanical energy to suppress cooling, it remains unclear 

According to observations (Jones & Forman M984I Fabian ... . ,. r, . . , b , , 

et al.[1994 Peres et al. 1998, Chen et al.lZOOTll. relaxed galaxy h ™ thlS ener ^ 18 dls ^ ted ln a homogenous way, such that 

, . . c . c a + u u-. . j. • 4. • -c cool cores appear in the form observed. Nevertheless, the CC 

clusters are most often found to exhibit at their centers a sigmfi- .^ XT • , , , 

. , ... T ™. . . , . .. ,. T . - AGN connection has been unambiguously demonstrated with 

cant drop in the ICM gas temperature due to radiative cooling. It . „ , . 5„„, ^ , , , 

• l . ii • i ,°, , iL . , ■ observations, e.g., Burns 1990, Eilek 2003, Sanderson etal. 2006 

is by now established that the cooling time in the central regions , f — , ' ; ,. .. — 1 

c , . • , . TT , , , ° • . . .. report that nearly every cool core cluster hosts a radio emitting 

of clusters is shorter than the Hubble time, thus originating very , fL, b 

, rp, c , c J* . AGN, creating cavities in the X-ray gas. 

dense cores. Therefore, we expect to observe a surface bright- ' " J b 

ness (SB) peak in these regions (with a typica l radiu s of the or- There is clear evidence for clusters harboring cool cores, 

der of 70 kpc - see for example Vikhlinin et al. 120031 along with however, many clusters do not show any signatures of cooling, 

some other features such as an enhanced iron abundance. At th is point, it may be pertinent to question whether there is 

The well known cooling flow problem stems from the ob- a CC/Non CC bimodal cluster population. In this the case, are 

servation that the few detected emission lines are not as strong CC clusters primordial and were some of them (the non-cool 

as expected to justify the predicted cooling rate. Further obser- core clusters) have been wiped out by mergers or an AGN heat- 

vational evidence is found in the ratio of the central tempera- ln 8 overshoot? Conversely, some authors support a scenario in 

ture with respect to the global cluster temperature (T central /T), whichnon-CC evolve to CC if no major mergers occur (O'Hara 

which remains at a factor -1/3 (Peterson et al. 120031 Bauer et al. et al. | 20061 >. 

The role of mergers in destroying cool cores has long been 

Send offprint requests to: J. Santos, jsantos@mpe.mpg.de debated (Fabian et al. 1984), both from observational results and 
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predictions from simulations. Currently, observations seem to fa- 
vor CC destruction through cluster mergers (Allen et al. 2001 
Sanderson et al. I2Q061 >. Simulations, however, yield ambivalent 
results. Using adaptive mesh refinement simulations, Burns et 
al. 120071 advocate that non-cool core clusters undergo major 
mergers which destroy embryonic cool cores while CC clusters 
grow at a slower rate without suffering early substantial mergers. 
However, an analysis of smoothed particle hydrodynamics sim- 
ulations of X-ray cluster mergers by Poole et al. 2006 , suggests 
that the heating of the ICM arising from mergers is not sufficient 
to prevent cooling. They argue that, in the KCDM scenario, the 
merger rate is too small to account for the local abundance of 
non-CC clusters whose disruptions occurs on short timescales. 
We note, however, that a good agreement in the overall prop- 
erties of observed and simulated CCs/NCCs is found by Burns 
et al., who were able to qualitatively reproduce the temperature 
structure observed in cool cores. 

The present observational census on the abundance of cool 
core clusters sets a fraction of Cool Core Clusters (CCC) in 
the local Universe in the range 50 -70 % (Peres et al. 1998 
Chen et al. 120071 1. As we move to increasing redshifts, re- 
sults on this topic become more sparse. Using spatially resolved 
spectroscopy, Bauer et al. 2005 presented a comparative study 
of cool cores at low and intermediate redshift (mean z=0.06 
and 0.22, respectively). The cool cores at intermediate redshift 
show the same temperature decrement, T/T cenlra i ~ 3-4, as the 
nearby CC's, and have a frequency rate similar to the local one. 
The central cooling time measurements confirmed this result. 
Consequently, they find no signs of evolution of the cool core 
fraction in the redshift range z ~ 0.15 - 0.4 to now. 

At present, the knowledge of the high redshift cluster pop- 
ulation is still rather poor: only 7 X-ray selected clusters with 
Z > 1 have been confirmed (Rosati et al. 1999, Perlman et al. 
120021 Stanford et al. [20021 Mullis et al. |2005l Stanford et al. 
2006). In particular, determining and quantifying the existence 
of cool cores at high redshift remains largely unexplored due 
to the lack of photon collecting power of current instruments to 
provide spectra with sufficient statistics for distant CC's. A first 
attempt to characterize a peaked core in high redshift systems 
was undertaken in Ettori et al. 2004, using double beta model fit- 
ting. The outcome of this analysis states that a second beta model 
component cannot be required to describe the distant cluster data 
- at least, not in the same way as observed for nearby cool core 
clusters. Recently, Vikhlinin et al. 2006 presented a study on the 
evolution of CCC at z > 0.5 based on pure imaging data, using 
a distant cluster sample drawn from the 400SD survey imaged 
with Chandra. These authors claim a strong evolution of the CC 
fraction in this redshift range, implying that we are missing cool- 
ing flows in the distant cluster population. However, this claim 
has been recently challenged with simulations by Burns et al. 
12007 who found no evidence for such a decline in the Cool Core 
fraction in simulated data with redshift up to one. 

The understanding of the formation and evolution of cool- 
ing flows is fundamental in pursuing cosmological studies using 
galaxy clusters. It is already known that cool core clusters affect 
the scaling relations causing a departure from the theoretical ex- 
pectations and increasing the scatter (Fabian et al. 119941 Zhang 
et al. 120061 ) . Whether or not we can put together a consistent pic- 
ture of clusters with and without cool cores will greatly affect the 
reliability of cluster results in deriving cosmological parameters. 

In this paper we investigate the fraction of cool core clusters 
in a high redshift sample imaged with Chandra at 0.7 < z < 1.4. 
To this aim we define and quantify criteria to select clusters as 
CCC in a sample of nearby clusters, 0.15 < z < 0.3, which 



are then applied to the distant sample. The low-z clusters exhibit 
markedly different morphologies which is essential to charac- 
terize cool cores with varying strength. The low photon statis- 
tics of the distant clusters made a spatially resolved tempera- 
ture analysis inviable, constraining our study to a spatial analy- 
sis of the cluster images. We analyzed the overall surface bright- 
ness and determined cooling times in the clusters' central re- 
gions. The surface brightness diagnostics were developed based 
on the nearby sample and were validated to high-z using sim- 
ulated distant clusters. These were obtained with a novel sim- 
ulation method which we present here, the cloning technique, 
in which we simulate distant galaxy clusters by emulating the 
appearance of the low-z clusters at the same redshifts of the ob- 
served distant sample. Such simulations allow us to make an un- 
biased comparison of cluster spatial properties between nearby 
and distant populations. 

The two cluster samples are drawn from the Chandra archive. 
The high-z sample, which is the most relevant for this study, 
is a statistically complete sample based on the Rosat Distant 
Cluster Sample (RDCS), from Rosati et al. [19981 This dataset 
allows us to draw conclusions on the abundance of CCC when 
the Universe was ~ 1/3 of its age. On the other hand, the nearby 
sample was constructed by taking morphologically interesting 
archive clusters and was used as a test set to understand how we 
can characterize cool cores. 

We note that archive samples will likely be biased in the 
sense that the majority of the observed clusters are "interesting", 
where interesting here means showing signs of mergers / disrup- 
tions, or on the contrary, having very regular shapes. In addi- 
tion, cluster selection based on flux limited samples will prefer- 
entially find high luminosity/mass clusters at high redshift. This 
effect, goes in the direction of favoring CC clusters (which have 
a bright core) with increasing z. For these reasons, constructing 
complete, unbiased samples to address issues of cluster evolu- 
tion presents difficulties which should be carefully dealt with. 
Since the distant and nearby samples are not drawn from a ho- 
mogenous sample with a known selection function, we do not 
address the issue of evolution of cool cores on this paper. 

As a complement to this X-ray study and with the aim of un- 
derstanding sample selection effects, we shortly investigate the 
surface brightness properties of an optically selected sample of 
high-z galaxy clusters imaged with Chandra. We assess the inci- 
dence rate of cool cores and compare it with the results on the 
X-ray selected cluster set. 

The paper is organized in the following manner: the cluster 
samples are introduced in section 2 and in section 3 the surface 
brightness profiles are presented, together with the beta model 
fitting. In section 4 we introduce the surface brightness concen- 
tration parameter and present the simulation method we used to 
test its application. In section 5 we present the cooling time anal- 
ysis. Section 7 is dedicated to a comparison of X-ray surface 
brightness properties in optically and X-ray selected high red- 
shift galaxy clusters. Finally, in section 8 we discuss the results 
obtained and conclude. 

The cosmological parameters used throughout the paper are: 
#0=70 km/s/Mpc, Q A =0.7 and Q m =0.3. 

2. Cluster samples 

We used archive Chandra ACIS-I and ACIS-S images in the [0.5- 
5.0] keV energy range: the 11 low-z clusters (Table 1) were se- 
lected from the sample compiled by Baldi et al. 2007, and the 
high redshift X-ray selected sample, which comprises 15 clus- 
ters (first part of Table 2), was drawn from the sample presented 
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Fig. 1. Nearby cluster sample observed with Chandra (point sources are removed). Images are smoothed with a gaussian of sigma=3 
and rescaled by the square root of the intensity. Individual boxes have a size of 8x8 arcmin. Cluster redshift increases from top-left 
to bottom-right. 
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Fig. 2. X-ray selected, distant cluster sample observed with Chandra (point sources are removed). Images are smoothed with a 
gaussian of sigma=3 and rescaled by the square root of the intensity. Individual boxes have a size of 3x3 arcmin. Cluster redshift 
increases from top-left to bottom-right. 



in Balestra et al. 120071 We therefore refer the reader to these 
two papers for details in the data processing. In Table 2 we also 
present the optically selected distant sample which will be dis- 
cussed in section 7. 

Point sources were removed from all images. The gaps be- 
tween the detector chips were filled assuming Poisson statistics 
and interpolating with neighboring pixels (dividing images by 
exposure maps is not enough for gap removal because the im- 
ages are subject to noise propagation). The local background 
was determined in empty regions of the cluster images; when 



computing the background a Gaussian filter with a kernel of 1 
was applied to the images to avoid zero pixel values. The cluster 
number counts were determined using the growth curve analysis 
(Bohringer et al. 120001) : cluster counts are calculated in circular 
apertures until the background level is reached. The maximum 
value of the growth curve was taken as the nominal cluster num- 
ber of counts. For each cluster we determined the center-of-mass 
coordinates, and assigned this value to the cluster center. 
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Table 1. Low redshift cluster sample: cluster name (1), redshift (2), spectral temperature (3), Right Ascension (4), Declination (5) 
and net exposure time (6) 



Cluster ID* L > 




T (keV) <J) 


RA (4 > 


DEC (S) 


Exposure (ksec) (b) 


A1413 a 


0.143 


7 cn+u.iO 
/,J -0.12 


11:55:18.20 


+23:24:28.80 


136 


A907 a 


0.153 


5.82 ±0.12 


09:58:22.00 


-11:03:50.40 


100 


A2104 a 


0.155 


6.76 ±0.19 


15:40:07.90 


-03:18:04.00 


48 


A2218 a 


0.177 


6.25 ±0.31 


16:35:56.00 


+66:12:45.00 


57 


A2163 b 


0.203 




16:15:45.77 


-06:08:55.00 


80 


A963 a 


0.206 


6 02 +a28 

-QiS 

7 4.^+0-49 

' -0.27 


10:17:03.40 


+39:02:51.00 


36 


A2261" 


0.224 


17:22:27.20 


+32:07:58.00 


32 


A2390 a 


0.228 


9.35 ±0.15 


21:53:36.50 


+ 17:41:45.00 


111 


A1835 11 


0.253 


8.09 ±0.53 


14:01:02.00 


+02:52:40.00 


30 


ZwC13146 a 


0.291 


8.59 ±0.39 


10:23:39.60 


+04:11:10.00 


45 


E0657 c 


0.300 


14.1 ±0.2 


6:58:28.60 


-55:56:36.03 


25 



M Temperatures taken from Baldi et al. 2007, [A) temperature from Markevitch et al. 120011 M temperature from Markevitch et al. 2005 

Table 2. High redshift cluster samples: cluster name (1), redshift (2), spectral temperature (3), Right Ascension (4), Declination (5) 
and net exposure time (6) 



Cluster ID ( 11 


z<*> 


T (keV) (j) 




RA' 4 ' 




DEC 5 ' 


Exposure (ksec) lb, 


RX J1221+4918 


0.700 


8 4 +u i! 


12 


21:24 


50 


+49:18:14.40 


78 


RX Jl 113-2615 


0.730 


5 6 +M 


11 


13:00 


00 


-26:15:49.00 


103 


RX J2302+0844 


0.734 


80^ 


23 


02:48 


06 


+08:43:54.72 


108 


MS 1137+6625 


0.782 


11 


40:22 


81 


+66:08:14.50 


117 


RX J1317+2911 


0.805 




13 


17:21 


84 


+29:11:17.00 


111 


RX J 1350+6007 


0.810 


4.5+ 07 


13 


50:47 


78 


+60:07:13.32 


58 


RX J17 16+6708 


0.813 


6 9 + ^ 


17 


16:48 


76 


+67:08:25.81 


51 


MS 1054-0321 


0.832 


75^ 

-°i 4 „ 


10 


56:58 


00 


-03:37:37.30 


80 


RX J 1226+3333 


0.890 


12.9+} | 


12 


26:58 


20 


+33:32:48.00 


9.5 + 31.5 


CL J1415+3612 


1.030 


7 +0 8 
/.o 07 


14 


15:11 


20 


+36:12:04.00 


89 


RX J09 10+5422 


1.106 


6.4!!- 


09 


10:45 


41 


+54:22:05.00 


170 


RX J 1252-2927 


1.235 


7 3 +l3 


12 


52:54 


50 


-29:27:18.00 


188 


RX J0849+4452 


1.261 


5 


08 


48:58 


52 


+44:51:55.08 


185 


RX J0848+4453 


1.273 


2 4 +25 


08 


48:36 


20 


+44:53:47.17 


185 


XMMU J2235-2557 a 


1.393 


6 +2 5 


22 


35:20 


70 


-25:57:40.70 


189 


RCS1419+5326 


0.620 


5 +U4 


14 


19:12 


14 


+53:26:11.56 


56 


RCS 1107.3-0523 


0.735 


4 3 +a5 


11 


07:24 


08 


-05:23:23.19 


93 


RCS1325+2858 


0.750 




13 


26:31 


04 


+29:03:25.02 


62 


RCS0224-0002 


0.778 


J -0.8 


02 


24:33 


61 


-00:02:24.68 


101 


RCS2318.5+0034 


0.780 


'• J -i.o 


23 


18:30 


67 


+00:34:03.03 


50 


RCS 1620+2929 


0.870 


4 6 +2 1 
5.3+J* 


16 


20:10 


01 


+29:29:15.41 


34 


RCS23 19.9+0038 


0.900 


23 


19:53 


59 


+00:38:09.05 


62 


RCS0439.6-2905 


0.960 




04 


39:37 


76 


-29:04:49.40 


92 


RCS1417+5305 


0.968 


14 


17:02 


13 


+53:05:23.57 


62 


RCS2 156.7-0448 


1.080 




21 


56:41 


63 


-04:47:53.47 


71 


RCS21 12.3-6326 


1.099 




21 


12:20 


51 


-63:26:13.97 


68 



M The temperature of XMMU J2235-2557 was taken from Mullis et al. 2005 Temperatures of RCS clusters refer to Bignamini et al. submitted. 
All other temperatures refer to the values given by Balestra et al. 120071 



2.1. Low-z sample 

The nearby cluster sample (0.15 < z < 0.3, median z=0.21) 
is presented in Fig. Q] The selection criteria for its construc- 
tion were essentially of observational nature since we were in- 
terested in low redshift clusters with a size that would properly 
fit within the field of view of Chandra (16 arcmin 2 ), and with 
sufficiently high number of counts (> 10 000) to ensure a proper 
analysis with good enough statistics of the cloned data (see sec- 
tion 4.2 for details on the cloning technique). The sample is het- 
erogeneous in the sense that we wanted to have clusters with 
known classification in terms of CC/non-CC; we therefore gath- 



ered this information from the literature - see Vikhlinin et al. 
120031 Bauer et al. [20051 R eese et al. 120021 De Filippis et al. 
120051 Arabadjis et al. 120021 Pratt & Arnaud 120021 According to 
these various sources, clusters A2104, A2218, A2163 and E0657 
are obvious non-coolcores whereas A1835 and ZWCL3146 are 
rare, high luminosity systems, exhibiting pronounced cool cores. 
Cluster temperatures obtained with spectroscopic fitting were 
taken mostly from Baldi et al. 120071 and also from Markevitch 
et al. |20TJTl & [20051 for clusters A2163 and E0657, respectively. 
Cores were excised when a cool core was present. All systems 
have T > 5 keV. 
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2.2. High-z sample 

The distant cluster sample (Fig. [2]i is drawn, mostly, from the 
RDCS (Rosati et al. I1998K with one important addition: a 
190 ksec Chandra exposure of XMMUJ2235. 3-2557 (hereafter 
XMMUJ2235), with z=1.393. This massive object is the high- 
est redshift cluster in our sample and the second most distant 
cluster (Mullis et al. 120051 1 confirmed to date (the most distant 
cluster known so far has redshift 1 .45 (Stanford et al. 2006 ). The 
sample ranges from 0.70 to 1.393 in redshift (median z=0.83), 
which spans a lookback time between 6.3 and 9.0 Gyrs, for 
the assumed cosmology. These are all massive clusters, span- 
ning approximately a decade in mass, with the exception of 
RXJ0848+4453 which has a temperature of 2.4 keV. The re- 
maining clusters in the sample have bulk temperatures which 
range from 4.4 to 12.9 keV and were obtained by spectroscopic 
fitting. Temperature measurements and general Chandra imag- 
ing reduction was extensively described in Balestra et al. 2007 
Concerning XMMUJ2235, we quote the temperature obtained 
by Mullis et al. 2005 using XMM-Newton data, although a more 
precise spectroscopic analysis of the new Chandra data will be 
available soon (Rosati et al. in preparation). Cluster number 
counts range from 250 to 11000 with a median value of 1750 
counts. 



3. Surface Brightness profiles 

3. 1 . Scaled Surface Brightness profiles 

Galaxy clusters are expected to exhibit self-similarity, which al- 
lows us to link their physical properties with temperature. This 
prediction is based solely on gravitational arguments where clus- 
ters are described as dark matter dominated entities. However, 
there has been evidence for a departure from the self-similar 
scaling scenario (Arnaud et al. 2005b, Chen et al. 2007): the 
Lx oc T 2 scaling relation expected in the case of gravitational 
only processes, is at odds with the observed L x oc r 2 88 rela- 
tion. This implies the intervention of non gravitational phenom- 
ena that provide extra energy input (Arnaud & Evrard 1999). An 
empirical scaling law was then derived which fits to the observa- 
tions. 

For the comparison of the local and distant clusters, we 
scaled the surface brightness profiles according to, both the em- 
pirical and self-similar scaling relations. In a self-similar, purely 
gravitational scenario, the surface brightness scales as Sx K 
r 05 , whereas in an empirical derivation Sx T lM (Neumann 
et al. [19991 Arnaud et al . 120021 . 

We account for the surface brightness dimming and consid- 
ered the redshift evolution of the Lx-T and R — T relations (see 
e.g. Ettori et al. |2004| l obtaining: 

SxKz)- 3 KT a (l+z) 4 K corr (1) 

where h(z), the cosmological factor, is defined as the Hubble 
constant normalized to its the current value, H z /Hq = 

■^Q. m (l +z) 3 + Oa- The slope a is dependent on the used scaling 
relation and K corr is the K-correction, which is defined in section 
4.2.1. 

The way in which scaling relations affect our analysis needs 
some clarification. The purely self-similar gravitational scenario 
(eg. Navarro, Frenk & White 1996, 1997), which seems to be 
reasonably well supported by observations, predicts an evolu- 
tion of the dark matter halos becoming more compact for fixed 
enclosed mass with increasing redshift. This is reflected in the 



global appearance of the ICM and the cluster morphology in X- 
rays, including further modification due to changes in the ratio 
of gas to total mass (Arnaud et al. 2002). This intrinsic evolution 
of the global cluster appearance is accounted for by the above 
mentioned empirical scaling relations. Focusing on the study of 
cool cores in clusters we are rather more interested in the local 
physical conditions at the cluster center which are for example, 
characterized by the cooling time of the gas or the total radiative 
energy loss of the ICM in the cool core region. This character- 
ization is not subject to the above described scaling relations. 
Thus if we are interested in comparing the rate at which the cen- 
tral ICM is radiatively cooling in nearby and distant clusters, we 
have to compare the cluster properties in fixed physical units. We 
will therefore in the following use the empirical scaling relations 
to reveal changes of the global cluster and cluster ICM structure 
properties, while unsealed relations are used to directly compare 
physical conditions in the cooling cores. 

We adopted a fixed density contrast A=500 with respect to 
the critical density of the Universe at the cluster redshift (p c (z)), 
to study the global properties of clusters. The fiducial radius, 
/?50o, was determined using the R^q - T relation in Evrard, 
Metzler & Navarro 1996, represented by the following fitting 
function: 

# 5 oo = h(zr l 2A8(kT/l0keVf 5 H 50 /H 70 (2) 

Eq. (2) is applicable for clusters with T > 3.5 keV. 

From the top panels in Fig. [3] it is clear that the low-z pro- 
files converge well for r > 0.2.Rsqq. However, at high-z two out- 
liers are evident in the self-similar scaling (bottom-right plot): 
MS1054-0321 (hereafter MS1054) and RX J1226+3333. The 
former is a merging cluster and the latter is the hottest distant 
cluster. Interestingly, RX J1226+3333 does not appear as an out- 
lier when using the empirical scaling. This just reinforces the 
better use of the empirical relation which does not boost the clus- 
ters at the extreme (highest temperature) outside the bulk of the 
data points in the relation. 

A striking difference between the CC and non-CC clusters 
in the nearby sample is evident in Fig. [3] (top panels). The flat- 
ness of the profiles of non-cool core clusters A2104, E0657, 
A2218 and A2163 are in obvious contrast with the remaining 
profiles which have a much higher central surface brightness - 
up to nearly 2 orders of magnitude larger - and a steepness which 
varies according to the strength of the cool core. This result is in- 
dependent of the choice of scaling. At high redshift this marked 
distinction is not apparent. 

The scatter of the surface brightness profiles scaled with 
T 5 amounts to 50% of the average value in the low-z sam- 
ple at r=0.3i?soo- In the distant sample the scatter is boosted by 
MS 1054; if we neglect this outlier we obtain a dispersion around 
the mean value of 65% at the same radius. 

3.2. Single and Double Beta model fitting 

Cluster surface brightness profiles are commonly analyzed with 
the approximation of the single isothermal beta model (Cavaliere 
& Fusco-Femiano 1976). We computed azimuthally averaged 
surface brightness profiles in bins of fixed radius, which were 
fitted with a single /3 model: 

S(r)=S„(l +(r/r c ) 2 r 3l3+0 - 5 + C (3) 

where So, r c , (3 and C are the central surface brightness, core 
radius, slope and constant background, respectively. The fitting 
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Fig. 3. Left Empirical scaling of SB profiles - low-z sample (Top) and high-z sample (Bottom) Right Self-similar scaling of SB 
profiles - low-z sample (Top) and high-z sample (Bottom). The red dashed lines corresponds to the averaged profiles. 



procedure was performed with a Levenberg-Marquardt least- 
squares minimization. 

When a cluster harbors a cool core, a single beta model is 
often inappropriate to describe its central excess emission, re- 
quiring the use of a 2B model: 

S(r) = Soid + (-) 2 )- 3A+0 - 5 + Soad + (-) 2 r %+0 - 5 + C (4) 

Such a conclusion has been remarked observationally (Jones 
& Forman 119841 Vikhlinin et al. 2005 ) and is also expected 
from simulations (Burns et al. 120071) . However, this observation 
cannot be so clearly stated for the high-z clusters because the 
low statistics allows to obtain a good x 2 re 4 using either a sin- 
gle or a double B function. This has been previously noted by 
Ota & Mitsuda 2004 who analyzed a sample of 79 clusters with 
0.1 < z < 0.8, and in Ettori et al. 120041 where nearly all dis- 
tant clusters in our sample were analyzed. There are however 
exceptions such as cluster CL J1415+3612 at z=1.03, in which 
the single beta model fails to capture the central SB. A double B 
model is qualitatively more appropriate to describe the core SB 
excess, although there are no significant differences in the re- 
duced^ 2 values of the IB and 2B model fit: ^ 2 ci/ = 1 .043 in the IB 
model fit and/j^^O.954 in the 2B model fit. These results make 
it obvious that fitting single/double B models to quantify the cen- 
tral surface brightness is not conclusive for the high-z sample. 
The /3-model parameters will not be used further in our study to 
characterize CCs, as we found them less useful than the follow- 
ing approaches. We nevertheless report our results in Table 3 for 



comparison with other work. Since the double-/? model results 
give, in general, a good representation of the data within the er- 
ror bars, the model fits also allow for an approximate check and 
reproduction of our analysis. 



4. Surface Brightness Concentration 

4.1. Introducing c SB 

In principle, the most stringent proof for the detection of a cool 
core is given by the temperature decrease in the core with respect 
to the bulk of the cluster. Unfortunately, the poor photon statis- 
tics does not permit a spatially resolved spectroscopic analysis 
of the high-z sample. 

Since a central surface brightness excess is a primary indica- 
tor (Fabian et al. 1 1 984t of the presence of a cool core, we eval- 
uate the core surface brightness in nearby clusters by defining a 
concentration parameter as the ratio of the peak over the ambi- 
ent SB. The aim of this approach is to use a single parameter to 
make a practical initial classification of cool cores. The apertures 
corresponding to the peak and the bulk were chosen to yield the 
largest separation between the CC and non-CC domains. The 
csb parameter was optimized using the low-z sample, varying 
the radius of the central peak and external radius. The optimal 
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Table 3. Single and double beta model fitting results: core radius, r c , slope B and reduced^ 2 . Errors are not presented when the 
parameter value is at the limit imposed by the fitting procedure. 



\j3 fit 2/3 fit 



cluster 


P 


r r (kpc) 


Xred 


A 




r cI (kpc) 


Pi 




r c2 (kpc) 


Xred 


A1413 


0.473+0.002 


56+1 


1.821 


0.755: 


t0.057 


70±5 


0.900 




286±2 


1.153 


A907 


0.515+0.002 


42+1 


6.360 


0.635: 


t0.152 


29±6 


0.706:1 


:0.045 


131+16 


1.439 


A2104 


0.566+0.014 


161+6 


1.335 






- 


0.542:1 


:0.018 


162±7 


1.447 


A2218 


0.841+0.150 


222+30 


1.220 






- 


0.606d 


:0.007 


159±4 


1.326 


A2163 


1.003+0.056 


349+15 


1.277 






- 


0.517:1 


:0.003 


184±2 


2.909 


A963 


0.516+0.003 


68+1 


1.165 


0.702+0.063 


76±7 


0.900 




301±8 


1.133 


A2261 


0.540+0.008 


66+3 


1.147 


1.200 




49+11 


0.559:1 


:0.006 


80±5 


1.090 


A2390 


0.474+0.001 


46+1 


17.037 


0.619 


±0.020 


49±2 


0.900 




411±5 


2.561 


A1835 


0.528+0.001 


29+0.3 


3.530 


0.699 


±0.050 


34±3 


0.719:1 


:0.016 


179+12 


1.803 


ZwC13146 


0.573+0.001 


2+1 


5.012 


0.639: 


t 0.156 


18±5 


0.667:1 


:0.023 


93+9 


1.224 


E0657 


1.000 


477+4 


1.864 






- 


0.900 




433+3 


1.919 


RX J1221+4918 


0.702+0.067 


262+31 


1.285 






- 


0.763d 


:0.062 


281+28 


1.321 


RX J1113-2615 


0.668+0.066 


107+17 


1.178 


1.000 




98±46 


0.800 




188+84 


1.048 


RX J2302+0844 


0.451+0.041 


73+16 


1.039 


0.470: 


t0.041 


68±8 


- 




- 


1.155 


MS 1137+6625 


0.807+0.041 


157+10 


1.152 


1.000 




11+25 


0.800 




55+14 


0.961 


RX J1317+2911 


0.537+0.094 


67+29 


1.039 


1.000 




52±38 


0.800 




175+53 


1.020 


RX J1350+6007 


0.561+0.095 


162+45 


1.085 


0.683 


±0.523 


36±27 


0.800 




341+48 


0.909 


RX J1716+6708 


0.740+0.062 


165+20 


0.988 


1.000 




28±6 


0.695:1 


:0.071 


154+25 


0.930 


MS 1054-0321 


1.000 


443+5 


2.826 








0.800 




360+5 


2.924 


RX J1226+3333 


0.710+0.053 


131+16 


0.851 


1.000 




23±11 


0.727:1 


:0.063 


136+21 


0.857 


CL J1415+3612 


0.552+0.026 


61+8 


1.043 


0.704: 


t0.094 


140±40 


0.800 




35+12 


0.954 


RX J0910+5422 


0.997+0.320 


177+57 


1.017 


1.000 




39+13 


0.800 




158+19 


1.022 


RX J1252-2927 


0.557+0.067 


89+21 


1.103 


1.000 




33+18 


0.601:1 


:0.111 


108+41 


1.114 


RX J0849+4452 


0.537+0.05 


72+19 


1.143 


1.000 




29±14 


0.800 




127+17 


1.043 


RX J0848+4453 


0.590+0.179 


122+69 


0.823 


1.000 




32±54 


0.800 




188+62 


1.193 


XMMU J2235-2557 


0.595+0.028 


94+11 


1.197 


1.000 




33+18 


0.800 




155+11 


1.090 



cs b is found for a peak radius of 40 kpc and a cluster bulk radius 
of400kpc: 

_ SB(r<A0kpc) 
° SB ~ SB(r<400kpc) ( ) 

To support this result we carried out extensive Monte Carlo 
simulations of CC and non-CC clusters modelled with, respec- 
tively, double and single beta profiles to optimize the choice of 
the internal and external apertures in the definition of c$b- We ex- 
plored the range of 30-80 kpc around the typical observed size of 
cool cores in nearby clusters. The chosen inner radius of 40 kpc 
is the one which make the separation between CC and non-CC 
clusters more obvious. The concentration parameter, as defined 
in eq. (5), was measured in these simulated clusters and we find 
that it efficiently separates the two classes of objects. 

Defining cs b in terms of fractions of a scaled radius such as 
^500 proved not to be suitable because the wide range covered in 
redshift implies sampling cluster regions of various sizes for dif- 
ferent redshifts. For example, 0. l/?5oo amounts to 1 32 kpc for the 
most nearby cluster, whereas for the most distant cluster it cor- 
responds to 65 kpc. Keeping in mind that the cool core region 
has been broadly described as the inner 70 kpc of a cluster, it is 
evident that using fractions of Z?5oo will introduce a bias towards 
lower redshift systems appearing more concentrated. In addition, 
scaling the radius with temperature would only be meaningful if 
the cluster samples spanned a different temperature range and/or 
if there were poor (T < 3.5 keV) systems, which is not the case 
here. Therefore we used a physical radius for measuring the con- 
centration index. 

Having thus defined the empirical cs b parameter using the 
low-z sample, we investigate whether there is a redshift depen- 
dence bias associated with this parameter which would inval- 
idate the comparison of c$b results between the samples. To 



this end we developed a simulation method which clones the 
observed low-z clusters to high redshift. Measuring the c$ b of 
simulated distant clusters provides a straightforward way to test 
for a redshift dependence. We describe this technique in the fol- 
lowing section. 



4.2. Testing c SB : simulating distant galaxy clusters with the 
Cloning method 

The cloning technique has been applied in an optical study by 
Bouwens et al. 1998, with the purpose of quantifying the evolu- 
tion of faint galaxies, using low-z galaxies in the Hubble Deep 
Field. By applying a pixel-by-pixel K-correction map to high 
resolution images of bright galaxies, a set of no-evolution deep 
fields was created, taking into account the space density and the 
cosmological volume. 

We have revised and adapted the cloning method with the 
purpose of simulating distant clusters of galaxies which are 
taken as standard non-evolving clusters that can be compared 
with distant samples. We clone the low-z sample to the redshifts 
matching those of the distant sample and for consistency in the 
comparison we normalize the simulated cluster counts to the cor- 
responding expected high-z cluster counts, adding as well the 
observed high-z background. Chandra angular resolution is es- 
sential for our analysis. Its point spread function (PSF) with a 
FWHM of ~ 1 arcsec corresponds to 8 kpc at z=\. Considering 
that the typical physical size for the CC region is about an order 
of magnitude larger than the instrument resolution, there is no 
need to deconvolve the images with the PSF which would make 
the application of our method very difficult. One of the main 
advantages of this technique is that it is model independent, im- 
plying that the simulations are under control, as no tuning of 
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parameters is required. The assessment of the selection effects is 
rather robust since we only have to rely on the data. 

The procedure comprises essentially 2 steps: the flux de- 
crease and the spatial rescaling which the template cluster un- 
dergoes. 

4.2.1 . Flux dimming: luminosity distance and K-correction 

The cluster flux decrease is due to both the physical luminosity 
distance, Di, and the redshifted energy band. In the adopted cos- 
mology Dl is given by (see for instance, Carroll, Press & Turner 
1992, Perlmutter et al. [T997l ): 



Ho Jo 



h(zy l dz 



(6) 




T=10 keV 



0.0 



0.5 



The X-ray flux is thus given by 



AnDj 



(7) 



where Lx is the X-ray luminosity in the observed frame. 
Therefore, when cloning to a higher redshift, we must compute 
the ratio of the squared luminosity distances. 

As we are using images in a given energy band (0.5-5.0keV), 
a bandpass correction - K-correction - plays an important role in 
the flux rescale because we are dealing with high redshift values 
(Hogg et al. |2002| l. If f y is the unabsorbed spectral flux density, 
we can express the K-correction for an energy band [E\, E<i\ as: 



Mzu^^- 2 )ecf(^.,^, Zl ,N H ,Z) 
Mz2,^,& 2 )ecf(^,^,z 2 ,N H ,Z) 



(8) 



where ec f stands for energy conversion factor, a term which re- 
lates flux to counts, and depends on the cluster redshift, hydro- 
gen column density, Nh, and metallicity, Z. 

Since the emissivity is not very temperature sensitive in the 
energy band used for the surface brightness measurements, it 
is practical and justified to treat the cluster ICM as isothermal. 
The K-correction terms are dependent mostly on the cluster red- 
shift (this dependence is stronger for z > 0.5) and temperature. 
The latter has a steeper effect for poor systems (T < 3 keV) 
which is not the case here. Applying single cluster temperatures, 
K-corrections were computed with XSPEC vll.3.1 (Arnaud ?). 
We used the optically thin MEKAL model spectra (Kaastra et 
al. [19921 with photoelectric absorption (WABS(MEKAL)), set- 
ting the metallicity to 0.3 solar and Af//=8. 10 20 cirT 2 . The spectra 
were folded with Chandra responses. The metallicity and the Nh 
have a weak influence in the K-correction: at the median red- 
shift and temperature of the high-z sample, using a higher value 
for the metallicity, e.g. Z/ZO=0.6, results in a decrease of the 
K-correction by x 8%, and a lower Nh value of 10 20 cirT 2 has 
the marginal effect of raising the K-correction by « 2% - see 
Fig- IH for a description of how the K-correction varies in these 
parameters. 

Eq. (8) thus becomes: 



[0.5-5.0]keV 



■£5 f* dv 
r 5.0(1 +z) 

Wa fvdv 



(9) 



Fig. 4. K-correction ratios computed with a MEKAL model: we 
show the dependence with redshift for temperatures [2,6,10]keV. 
Red lines refer to an Nh of 8.10 20 cm~ 2 and a metallicity of 0.3 
solar; blue lines refer to Nh = 8.10 20 cm' 2 and Z/Zo=0.6 and 
green lines refer to N H = 1.10 20 cm' 2 and Z/Z©=0.3. 



minimum of 1.5, when we clone a high temperature cluster at 
z ~ 0.3 to a final redshift of 0.7, and may reach a maximum of 4 
when the cloning redshift is the largest value, i.e, 1.393. 

If a second temperature component related to the cluster core 
were considered, this would increase the K-correction by a factor 
of 10-20%. 



4.2.2. Spatial rescaling 

The spatial rescale which a cluster undergoes when cloned to a 
higher redshift corresponds to the ratio of the angular distances, 
Da - Dl/(1 + z) 2 . Since the premise of this method is non- 
evolution of the clusters, it should be clear that the size rescale 
does not account for the shrinking expected when applying the 
cosmological factor h(z), which comes about as an intrinsic evo- 
lution of clusters. The spatial resizing was implemented using 
a gridding linear interpolation. Note that, as explained above, 
apart from the redshift resizing, no additional size scaling of the 
kind of eq. (2) was applied, as we compare the cluster cores in 
physical units. 

To illustrate this method we show in Fig. [5] an example of 
cloning. We clone A2163, a non-cool core cluster at z=0.2, to 
redshifts 0.7 and 1.0. The observed high-z clusters at the same 
redshifts (RXJ1221 at z=0.7 and CLJ1415 at z=1.0) are dis- 
played in the lower panels. The simulated X-ray images were 
normalized to have an identical number of counts as the cor- 
responding high-z distant clusters, and the same background, 
which were measured as described in section 2. 

A simple consistency check which allows us to evaluate the 
procedure resides in measuring the surface brightness within a 
physical aperture, both in the real and simulated image. For a 
given energy band we should obtain a ratio of the integrated flux 
which scales as follows: 



Sx(zi) = (1+z/) 4 
Sx(zf) (1+z,-) 4 



(10) 



so that Lx - Lxo/^io.5-5.0], where Lxo is the rest frame luminos- 
ity of the source. The calculated correction ratios range from a 



We recover perfect agreement of the results from the two 
approaches (within the expected numerical errors). 
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Clare, z=0.7 Clone, z=1 .0 




V1 221. 1=0.7 CIJ1415, z=1.0 



Fig. 5. Cloning example: Chandra image of low redshift non- 
CC cluster, A2163, and 2 cloned images (top), compared with a 
non-cool core cluster at z=0.7, RXJ1221 (bottom left), and a CC 
cluster at z=1.0, CLJ1415 (bottom right). Green circles have a 
radius of 2 arcmin. 

4.3. The c$b distribution 

We find a very good agreement between the surface brightness 
concentration of the clones and the parent (low-z) c$ b distribu- 
tion, implying no redshift dependence in c$b- This can be ex- 
emplified by making 10 realizations of cloning the complete 
nearby sample to 0.7 < z < 1.4 and computing the ratio 
csB[clones]/csB[low - z] as a function of redshift. This is illus- 
trated in Fig. [6] the ratio of c$b randomly fluctuates around one 
with no trends in redshift, which is confirmed by the weighted 
mean linear fit to the data. Moreover, no dependence with tem- 
perature was found. This indicates that cs b is an unbiased quan- 
tity which can be applied to the real high-z sample, therefore we 
can safely compare the nearby with the distant samples using 

CSB- 

The distribution of c$ b for both the nearby and distant clus- 
ter samples in presented in Fig. [7] the histogram of the low- 
z sample is shown in blue with the CCC in blue-hatched and 
the high-z sample is represented in red. A Kolmogorov-Smirnov 
(KS) test yields a probability of 45% that both csb distributions 
were drawn from the same population, implying there is no sig- 
nificant difference between the two samples with respect to the 
cs b parameter, cs b was measured by performing aperture pho- 
tometry in central regions enclosing 40 and 400 kpc. Errors were 
determined using Poisson statistics - see Table 4 for a listing of 
these values. 

5. Stacked SB profiles in bins of c S b 

Owing to the low number of counts of the distant cluster sample, 
we resorted to stacking the scaled SB profiles which we com- 
puted in section 3.1, making use of the csb results (Table 4). 
From Fig. [7] we define three categories of cool cores: non-CC 
(csb < 0.075), moderate (0.075 < c$b < 0.155) and pronounced 
(c SB > 0.155) CC. 

Fig.[8]shows the stacked profiles of both samples scaled ac- 
cording to the empirical scaling law, in units of Rsoo- At high-z 
we are limited by the spatial resolution of Chandra up to 0.03 
7?5oo, and at low-z we are constrained by the background and 
cluster size to probe regions up to 0.5i?soo (with the exception of 
A2104 which can be traced only as far as 0.3 Rsoo). The clusters 





[-t---= 
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redshift of clones 



Fig. 6. Analysis of 10 realizations of cloning the nearby clus- 
ter sample to high-z: we plot the ratio csB[clones]/csB[lowz] as 
a function of cloning redshift. Individual error bars refer to 1- 
sigma confidence level. The red solid line represents the linear 
fit to the data and corresponding 1-sigma errors are shown in 
dashed line. 



Table 4. Surface brightness concentration, c s B of the nearby and 
distant samples 



Low-z ID 




s 


High-z ID 




j 


A1413 


0.095 ± 


0.001 


RXJ1221+4918 


0.028 ± 


0.004 


A907 


0.169 ±0.002 


RXJ1 113-2615 


0.095 ± 


0.012 


A2104 


0.044 ± 


0.001 


RXJ2302+0844 


0.072 ± 


0.009 


A2218 


0.042 ±0.001 


MS 11 37+6625 


0.096 ± 


0.007 


A2163 


0.024 ± 


0.002 


RXJ1317+2911 


0.123 ± 


0.030 


A963 


0.101 ±0.002 


RXJ 1350+6007 


0.057 ± 


0.012 


A2261 


0.111 ± 


0.003 


RXJ1716+6708 


0.082 ± 


0.010 


A2390 


0.120 ± 


0.001 


MS 1054-0321 


0.016 ± 


0.002 


A1835 


0.236 ± 


0.003 


RXJ 1226+3333 


0.083 ± 


0.011 


ZwC13146 


0.217 ± 


0.003 


CLJ1415+3612 


0.151 ± 


0.015 


E0657 


0.017 ± 


0.001 


RXJ09 10+5422 


0.101 ± 


0.021 








RXJ1252-2927 


0.088 ± 


0.014 








RXJ0849+4452 


0.099 ± 


0.023 








RXJ0848+4453 


0.064 ± 


0.030 








XMMUJ2235-2557 


0.103 ± 


0.012 



with complex morphology due to mergers (E0657 at low-z and 
MS 1054 at high-z) were excluded from the stacking procedure 
because they introduce strong deviations from the average pro- 
file, which are related with the difficulty of identifying the cluster 
center. The stacked profiles exhibit different shapes (slope and 
central emission) consistent with Fig. [7] where the low-z sample 
spans over 3 bins of cs b whereas the high-z sample covers the 
first 2. Both the low and high redshift profiles are in full agree- 
ment on the outskirts at r « 0.3/?5oo (the deviation of the high-z 
NCC profile (binl) at large radius is mostly due to noise of the 
individual profiles), but not at the centers where we find a central 
emission offset: at r/R 500 = 0.03 the SB of binl is 0.023+0.005 
and 0.014+0.006 cts/s/arcmin 2 for the low-z and high-z samples 
respectively, with a ratio SB(low-z)/SB(high-z)=1.6. At the same 
radius the offset is larger for the second bin: SB=0. 089+0.026 
cts/s/arcmin 2 for nearby clusters and 0.029+0.019 cts/s/arcmin 2 
for distant clusters, with a corresponding ratio of 3.0. The quoted 
errors refer to the standard deviations associated with the stack- 
ing process, which are larger than the combined measurement 
errors. 

The central surface brightness clearly does not follow the ex- 
pected predictions from scaling relations. As mentioned above, 
it may be more meaningful to compare the physical state of the 
ICM directly, thus droping the temperature and redshift (h(z)) 
scaling. In this case, we obtain a closer match between the low- 
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Fig. 7. Histograms of c$ b for the low-z sample (blue solid line 
with the known CC clusters indicated in hatched) and for the 
high-z sample (red dash line). 
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Fig. 8. Stacked surface brightness profiles according to c$ b bins: 
the nearby sample (blue) presents 3 bins: non-cool core (solid, 
bin 1), moderate CC (dash, bin 2) and strong CC (dash-dot, bin 
3); the distant sample (red) shows 2 bins: non-cool core (solid, 
bin 1) and moderate CC (dash, bin 2). 

Z and high-z profiles, indicating that the physical conditions are 
similar. 

6. Cooling time analysis 

The radiative cooling of the intracluster gas originates mainly 
from thermal bremsstrahlung emission, with additional line 
and recombination radiation. The cooling time is defined as 
dt/dlnTgas (Sarazin 1988) and allows the evaluation of the cool- 
ing rate in galaxy clusters. Adopting a gas enthalpy model for 
the cooling function, t c00 i can be approximated by: 

2.5n g T 

tcool = 9 , ™, (11) 
njA{T) 

where A(T), n g , n e and T are the cooling function, gas num- 
ber density, electron number density and temperature respec- 
tively (Peterson & Fabian 120051 1. and with n g =\.9n e . The cen- 
tral cooling time is a sensitive parameter to characterize a cool 



core: when a cool core forms, the central temperature decreases 
and conversely the central density increases, resulting in a small 
cooling time. Thus we expect that CCC present cooling times 
shorter than the Hubble time or the time since the last major 
merger event. 

6.1. Density profiles 

Gas density profiles were obtained by deprojection of the sur- 
face brightness profiles. We corrected the gas density for the 
cosmological expansion: attending the M gas - T and R — T scal- 
ing relations (Arnaud et al. 2005a) it follows that n gas should be 
scaled by /z(z)~ 2 to account for the more compact shape of distant 
clusters. The density profiles of both cluster populations are pre- 
sented in Fig. [9] Similarly to what we found earlier in the scaled 
surface brightness profiles, the high-z density in the core region 
is also systematically lower than the central low-z gas density. 
This should be due to the h(z) scaling, as explained earlier in 
section 5. 




0.001 0.010 0.100 1.000 

r/r500 

Fig. 9. Scaled density profiles: low-z (blue) and high-z (red) clus- 
ters 



6.2. Cooling time at r=20 kpc 

The cooling time was measured within a central aperture of 20 
kpc radius extrapolating the density profiles with a 1/3 model. We 
found it preferable to use an aperture with fixed physical radius 
to make a consistent analysis of all clusters, as we are limited 
by the instrument's resolution at high redshift. Using a fraction 
of a scaled radius was not suitable as the redshift range is fairly 
large, and therefore we would probe larger areas at low redshift, 
introducing a bias. As a result, the strong cool cores would have 
a (artificially) boosted cooling time by a factor ~ 2-3. 

Since temperature profiles cannot be derived from high-z 
data, we used the spectral fit temperature in the single temper- 
ature (MEKAL) model for the cooling time computation. As a 
result, the t coo \ values in this case are an upper limit because 
in the presence of cool cores the central temperature would be 
lower, decreasing the cooling time. 

The cooling time distribution is plotted in Fig. [TOland pre- 
sented in Table 5; errors were obtained from the propagation 
of the errors on the temperature and ng. Four low-z clusters 
(A1835, ZwC13146, A2390 and A907) present a cooling time 
< 3 Gyrs and the 4 nearby non-CC show a t c00 i > 10 Gyrs. At 
high-z there is no system with t c00 i < 3 Gyrs although the major- 
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Table 5. Cooling time, t coo \ measured at 20 kpc, of the nearby 
and distant samples 



10" 



Low-z ID 



tcooi (Gyr) 



A1413 
A907 
A2104 
A2218 
A2163 
A963 
A2261 
A2390 
A1835 
ZwC13146 
E0657 



5 5$*^ 

J - Ju -0.17 

2.91 ±0.10 
15.10 ±0.99 
12.10 ± 1.22 

15.10!$ 
4.5231 
3.8S 

2.67 : 

1.13 ±0.09 

1.32 ±0.08 
20.80 ± 1.62 



High-z ID 



-0.26 
:0.07 



RX J1221+4918 

RX Jl 113-2615 
RX J2302+0844 

MS1137+6625 
RX J1317+2911 
RX J1350+6007 
RX J1716+6708 

MS 1054-0321 
RX J1226+3333 

CL J1415+3612 
RX J0910+5422 

RX J 1252-2927 
RX J0849+4452 
RX J0848+4453 
XMMU J2235-2557 



tcooi (Gyr) 
20 60 +4 35 
7 01 +1 

-OB 

14 70+ 3 Ss 

-3,67 

6.38^5 
8.77 + || 
10.90!|| 
6.39 + ™ 
16.70!|| 

4.11+0-88 

3 90 + 
8.65 + | i[ 
7.35 + ™ 
9.02+*jf 
6.19 +787 

4 84 + ^ 



ity (11/15) has t coo \ < 10 Gyrs, with only 4 clusters presenting a 
t coo i > 10 Gyrs. 

While similar cooling times or cooling rates in the nearby 
and distant samples point to similar physical conditions in the 
central cluster regions, the different ages of the clusters lead to 
some different conclusions of these results. The cooling rates 
for more than half of the nearby clusters are shorter than their 
estimated ages (~10 Gyr) which implies that they approach a 
steady state of cooling, some mass deposition, mass inflow and 
reheating by central AGN. With the same cooling time distri- 
bution, only a small fraction of the distant clusters (with ages 
between 5-8 Gyr) will have a chance to reach this steady state. 
If we therefore use the classical definition of cooling flows in 



fraction of CCs is 4/15, which is significantly smaller than the 
ratio found for nearby clusters (e.g. Peres et al. 1998, Bauer et 
al. 120031 Chen et al. 120071 ). 
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Fig. 10. Distribution of the cooling time, t coo i, in the low-z (blue) 
and high-z (red) samples. 



6.3. The correlation t coo i - c SB 

As we are interested in evaluating the degree of cooling of the 
central intracluster gas with a simple function that can be de- 
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Fig. 11. Correlation between cooling time, t coo i and the surface 
brightness concentration, csb, for the low-z (blue diamonds) 
and high-z (red squares) samples. The black solid line refers to 
the composite (both populations) bootstrap fit; the blue and red 
dashed lines refer to separate fits to the nearby and distant sam- 
ples, respectively. 



duced from imaging data only, we investigated the relation be- 
tween the cooling time and the surface brightness concentra- 
tion. We find a strong negative correlation with a Spearman's 
p rank correlation coefficient of -0.84 and a significance of non- 
correlation p - 8.15xl0~ 8 . Using 1000 bootstrap samples we 
performed a linear fit in the log-log plane of t co „i and c$ b, from 
which we derived an average slope of the power law function 
which describes the correlation between the two cooling esti- 
mators. A composite fit to both cluster samples yields t coo i oc 



which t coo i < toge also for the distant clusters, we find that the sb 



-1.10±0.15 



, (Fig.rTTIi. Fitting the nearby and distant samples sep- 



-0.92±0.26 
L SB 



arately yields the following correlations: t CO oi{high z ) • 
and tcooi(low z ) cc c^ 16±014 . 

Using the single fit correlation we obtain the cooling time at 
the reference thresholds c$b — 0.075 and csb — 0.155, as 6.9 
Gyrs and 3.1 Gyrs respectively. 



7. X-ray versus Optically selected distant clusters: 
a comparison using c SB 

In this section we briefly analyze the X-ray surface bright- 
ness properties of distant galaxy clusters detected in the optical. 
Specifically, we evaluate the degree of cooling using the previ- 
ously defined c$b parameter and compare with the results ob- 
tained for the X-ray selected high redshift sample. 

The Red Sequence Cluster Survey (RCS) aims at finding and 
characterizing distant galaxy clusters using R and z-band imag- 
ing, up to redshift z=1.4 (Gladders & Yee 2004). Selecting clus- 
ters in the optical allows one to study high redshift poor systems 
with low temperature, which are usually not found in flux limited 
samples restricted to high Lx objects at high redshift. 

Although there are indications of a large scatter in the Lx-T 
relation of the RCS sample, such clusters with low temperature 
(median T = 5 keV) enable us to probe a larger dynamical range 
in the Lx-T relation. 

The connection between optical and X-ray characteristics of 
these distant clusters has been under investigation (Hicks et al. 
120071 ). Currently, 11 clusters with 0.6 < z < 1.2 have been ob- 
served with Chandra and X-ray properties as temperature and lu- 
minosity have been derived (Hicks et al. 120071 Bignamini et al. 
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submitted). We based our analysis on Chandra images reduced 
by Bignamini et al. - see Table 2 for an overview of the sample. 

We find that 3 of the 1 1 RCS clusters are barely detected in 
X-rays and have such a low surface brightness that a c$b anal- 
ysis is not possible. The measurement of the surface brightness 
concentration as described in section 4.3, yields that the major- 
ity (6 of 8) of the RCS clusters with measurable X-ray emission 
have csb < 0.075, thus falling in the non-CC regime. Of the 
remaining two: RCS 1107 has a csb = 0.143 indicating a mod- 
erate CC and RCS1419 has c$b = 0.185, suggesting a strong 
CC - see histogram in Fig. [12] The c$b distribution of the X-ray 
distant clusters is overplotted in the same figure (red line) for 
direct comparison. The Kolmogorov-Smirnov statistic between 
the two cs b distributions is 0.4, with a probability of the samples 
being drawn from the same distribution of 20%. With such low 
statistics, the KS test does not provide a sensitive distinction on 
these two datasets. 

The outcome of this analysis based on our small samples 
suggests that high-z optically selected clusters may have a lower 
fraction of cool cores with respect to clusters detected in X-rays. 
Larger samples are needed to draw conclusions on possible dif- 
ferent physical conditions of the ICM in optically selected clus- 
ters, which may have not reached a final state of virialization 
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Fig. 12. Histogram of c$ b for the RCS sample in solid blue line; 
the high-z X-ray sample is overplotted in red dashed line. 

8. Discussion and Conclusions 

In this paper we investigate the detection and quantify the 
strength of cool cores in the high redshift cluster population. 

We obtained azimuthally averaged SB profiles and scaled 
them, in this manner accounting for different cluster tempera- 
tures. The segregation of cool cores at low redshift using this 
method is evident in Fig. [3] (top panels). However, this segre- 
gation is not obvious in the distant sample, and the use of sin- 
gle/double (3 model fitting does not provide conclusive results. 

A surface brightness concentration index, cs b, was then de- 
fined and measured in the local sample, which efficiently sep- 
arates the non-cool core from the cool core regime. This is in 
agreement with previously published results on the cuspiness of 
these clusters. In Hashimoto et al. 2007, a similar surface bright- 
ness concentration was defined, using the ratio of elliptical aper- 
ture radii enclosing, respectively, 30 and 100 % of the cluster 
surface brightness. This parameter was applied to a sample of 
101 clusters with 0.05 < z < 1, these including 4 of the 11 
nearby clusters and 12 of the 15 distant clusters we study here. 



Comparing the results from both measurements for individual 
clusters we find a good agreement with our results, with the ex- 
ception of A1835. 

We also developed a robust parameter free method to simu- 
late galaxy clusters at redshift z > 0.7, making use of the het- 
erogeneous nearby (0.15 < z < 0.3) sample. These simulated 
clusters serve as a benchmark for non evolving clusters, which 
allow us to test the suitability of indicators of cool cores at high 
redshift. Using our cloning technique, we verified the redshift 
independence of c$ b, thus making it a simple, unbiased quantity 
to study cool cores in the distant cluster sample (Fig. [7). The 
analysis of the concentration and radial profiles allows us to de- 
fine three csb bins which distinguish different regimes: non-CC 
(c SB < 0.75), moderate CC (0.75 < c SB < 1.55) and strong CC 
(csb > 1-55). By stacking the profiles in these bins (Fig. [8), we 
obtained a robust classification of cool cores which indicate the 
range of different cooling rates which the two cluster samples 
span. This is particularly useful for the distant sample which is 
affected by low X-ray count statistics. 

From the overall surface brightness analysis we conclude 
that the majority (10/15) of the high-z clusters present mild cool 
cores, similar to those found in the nearby sample. 

In addition to the surface brightness analysis, we measured 
the clusters' cooling time in a central region with 20 kpc radius, 
using a single global temperature. Analysing the reference, low- 
z sample, we conclude that strong cool cores present t coo i < 3 
Gyrs, whereas moderate CC have 3 < t coo \ < 10 Gyrs and finally 
non-CC show t coo \ > 10 Gyrs. The majority (11/15) of the high- 
z clusters are characterized by 3 < t coo \ < 10 Gyrs with only 4 
systems showing t coo \ > 10 Gyrs. We found no distant clusters 
with t coo i < 3 Gyrs. Similar cooling time bins were provided by 
Bauer et al. 2005 to discriminate between strong, mild and non- 
cool core systems. These authors measured the central cooling 
time and also the cooling time at 50 kpc: good agreement on their 
t coo \ at 50 kpc and our t coo i at 20 kpc is found for the common 
objects in the intermediate-z sample of Bauer et al. 2005] (4/6 
objects within lcr errors). Even so the distribution of the cooling 
time is quite similar in the nearby and distant samples, a smaller 
fraction of distant clusters has t coo i < t age ~ t Hubble, and therefore 
the fraction of distant CCs following the classical cooling flow 
distribution is smaller than for local clusters. 

We find a strong correlation between t coo i and c$ b, which al- 
low us to reliably relate a physical quantity with a phenomeno- 
logical parameter. This correlation is described by a power law 
fit: tcooi 00 Cgg *" (Fig. ITTb. The correlation does not change 
significantly if we fit the two samples separately. The low-z clus- 
ters show a slightly steeper slope (-1.16 + 0.14) when compared 
with the high-z slope (-0.92 + 0.26). 

Investigating galaxy cluster properties at high redshift is dif- 
ficult as these systems have small angular size, requiring high 
resolution instruments and long integration time, particularly for 
spectroscopy. It is therefore understandable that few attempts 
have been made to probe the centers of distant galaxy clusters. 
Vikhlinin et al. 2006 presented a study on the evolution of CCC 
in the redshift range [0.5-0.9] which is compared to a local sam- 
ple. They find evidence for evolution of the CC fraction, with 
a lack of cool cores in the distant sample. They use a different 
indicator for the surface brightness cuspiness which might have 
a lower sensitivity to moderate cool core systems. A meaningful 
comparison with our study would require a detailed analysis of 
the different cool core estimators. 

One often finds that massive, luminous galaxies containing 
AGN, lie at the centers of cool core clusters (Eilek 2003). In 
an attempt to establish such a correlation we searched for cD 
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galaxies in the cluster samples, examining available optical im- 
ages. We find that, in the low-z sample, all CCC indeed host a 
bright, massive galaxy, but also A2104, a non-cool core clus- 
ter, possesses a cD galaxy with an associated AGN (Liang et al. 
2000 Martini et al. 120041 1. Could this be a sign of an AGN heat- 
ing overshoot which could have contributed to destroy the cool 
core? At high redshift there is no obvious trend: the majority of 
the clusters possess a cD galaxy but it does not correlate with 

CSB- 

In the established scenario of hierarchical structure forma- 
tion, galaxy clusters develop through mergers and by accretion 
of neigbouhring smaller objects. Numerical simulations (Cohn 
2005 ) predict higher cluster merger rates with increasing red- 
shift, where the fraction of clusters with recent mass accretion 
due to mergers can be doubled at z=0.7, with respect to the local 
abundance. This framework provides a possible mechanism for 
preventing the formation of prominent cool cores at high red- 
shift. However, since high-z clusters are younger, other mecha- 
nisms may cause a delay in the formation of cool cores due to 
some internal energy release, such as AGN activity and star for- 
mation processes. The observed absence of pronounced CC at 
Z > 0.7 is therefore plausible in the current cosmological frame- 
work. 

We also investigated possible sample selection effects in the 
search for high redshift cool cores by studying 1 1 optically se- 
lected clusters from the RCS sample, which have a median tem- 
perature of 5 keV. We found that, unlike the X-ray sample, the 
majority of clusters in the RCS subsample have c$b < 0.075, i.e. 
they lie in the non-cool core category. 

To further extend this analysis we need to construct a com- 
plete low redshift sample, possibly with the same selection func- 
tion as the distant sample, to be able to trace the evolution of 
the abundance of cool cores. A good prospect for this study is 
the local representative cluster sample XMM-Large Programme 
(Bohringer et al. 120071 1. which comprises 33 clusters in the red- 
shift range z=0.055 to 0.183, drawn from the REFLEX sample 
(Bohringer et al. 2001). An intermediate redshift sample (z=0.3- 
0.6) from a similar project is also available (RI. M. Arnaud) 
which provides an intermediate step in the assessment of the cool 
core evolution with redshift. 

Clearly, having a high-z data set covering a larger volume 
would also be very advantageous. Serendipitous cluster surveys 
currently underway, such as the XCS (Romer et al. 1200 IK the 
XDCP (Mullis et al. [20051 Bohringer et al. [20051 and XMM- 
LSS (Pierre et al. [2006 ), will yield suitable distant samples for 
these studies in the near future. 

A comparison of our results with cosmological simulations 
would also provide additional new constraints to cluster forma- 
tion scenarios. In fact, new insights on the evolution of simu- 
lated CCs and NCCs have recently been published in Burns et 
al. 120071 which are in good agreement with our findings. 
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